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A quasi-particle model is employed to derive from available lattice QCD calculations an equation 
of state useable in hydrodynamical simulations of the expansion stage of strongly interacting 
matter created in ultra-relativistic heavy-ion collisions. Various lattice results give an astonishing 
agreement of the pressure as a function of energy density at large energy densities supposed 
the pseudo-critical temperature is in the range 170 ± 15 MeV, while in the transition region the 
equation of state is not yet well constrained. Therefore, one can construct a family of equations of 
state by bridging the uncertain region from the uniquely given high-energy density region part to 
a hadronic equation of state by suitable interpolation together with the extrapolation to non-zero 
baryon density by means of the quasi-particle model. We present a series of tests of the model, 
discuss the chiral extrapolation and the role of Landau damping. We also briefly sketch the path 
of cosmic matter in the early universe in the phase diagram. 
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1. Introduction 

The equation of state of strongly interacting matter is a central issue for understanding and 
modelling the adiabatic path of cosmic matter of the expanding universe (hot QCD), the expan- 
sion of matter created in ultra-relativistic heavy-ion collisions (hot and medium-dense QCD), and 
compact stars (dense QCD). Intimately related is the phase diagram of strongly interacting matter 
with challenges like the nature and localization of the deconfinement transition and the occurrence 
of a critical point on the phase border curve. Ab initio calculations evaluating observables which 
quantify these notions are still in progress. It is, therefore, opportune to employ at this stage of 
insight in QCD appropriate models to interpolate and extrapolate the various pieces of knowledge 
aiming at delivering, e.g., a useable equation of state. 

The tool employed here is the quasi -particle model developed in [[]]]. We are going to ad- 
dress the question whether a unique equation of state, applicable for ultra-relativistic heavy-ion 
collisions, is at our disposal. Furthermore, we report on a few recent improvements of the model. 

Our paper is organized as follows. In section 2 we review the quasi-particle model, which is 
exploited in section 3 to derive an equation of state suitable for predictions for heavy-ion collisions 
at RHIC and LHC energies. The possibility to extend the model by explicitly accounting for the 
critical point is addressed in section 4. Partial improvements (such as an attempt of chiral extrap- 
olation, and inclusion of Landau damping and other plasma excitations, and imaginary chemical 
potential) are briefly reported in section 5. Finally, we comment on the adiabatic path of strongly 
interacting matter in the expanding universe (section 6). 



2. Quasi-particle model 

The present basic version of the employed quasi-particle model is related to QCD as follows 
[§] : (i) two-loop <J> functional which results in one-loop self-energies, (ii) neglect of imaginary parts 
of self-energies (and, via Dyson's relation, also in propagators) as well as neglect of (anti)plasmino 
and longitudinal gluon excitations, (iii) use of approximate energy (ft)) - momentum (k) relations 
of quasi-particle excitations 



CO^=k 2 + mi, mi 



X - ([2N c +N f ]T 2 + ^N ff A G 2 (T,ii q ), (2.1) 



a$ = k 2 + m 2 q + 2m q M + + 2M 2 + , M 2 + = I T 2 + j G 2 (T, n q ) (2.2) 

for transverse gluons (T) and quarks (q), where (iv) the effective coupling G 2 obeys Peshier's 
equation a T ^- +a^ q ^- + acG 2 = with coefficients a^ TG given in [§, |], once G 2 (T,\x q = 0) 
in the parametrization 

G| oop (0, S=X^,T>T C 
G l {T) ={ (2.3) 

G|i oop (r c )+fc(i-|), t<t c 

is adjusted to lattice data. (Gfioop is the two-loop running QCD coupling, however, with argument 
(§.) The expressions for baryon and entropy density look like the standard statistical integrals, 
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however, with state-dependent mass gaps m 2 ^ and M 2 = m q + 2m q M + + 2M\ for gluons and quarks. 
For quarks, the existence of a mass gap related to M + emerges now from lattice simulations in the 
quenched approximation iQj. The dependence of the latter ones on the quark chemical potential 
jl q , beyond the dependence on the temperature T, is called the BKS effect in 

The model has been tested against various sets of lattice QCD results for Nf flavors and N c = 3 
colors, see [jj], g]. A recent extension of the model towards two independent chemical potentials 
\l u ^ allows the calculation of quark number and isovector susceptibilities as well as diagonal and 
off-diagonal susceptibilities and their respective Taylor expansion coefficients 

h=2c, + UcJ^\ 2 



P = 2c2 + 12c 4 ^J +•••, (2.4) 
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§=24 + 12d(|)%- ; (2.5) 



|f =2cr + 12cr(^) 2 + ---, (2.6) 



^- = 2cf + l2cT{^) +••■, (2.7) 

being second-order derivatives of the grand thermodynamical potential. The light-quark chemical 
potentials are decomposed as \i q = ^(fl u + fid) and = \{fl u — A t rf)> where jii denotes the isospin 
chemical potential. A comparison with lattice QCD results is displayed in Fig. [jj This additional 
test provides further confidence in the model. 

3. Hydrodynamics for RHIC and LHC 

Due to various differences of the implementation of QCD on a space-time grid to evaluate the 
equation of state for Nf ■ = 2 + 1 flavors, some differences are obvious, as demonstrated in the left 
panel of Fig. ||. The surprise however is that a translation of these differing lattice results for the 
scaled pressure p as a function of the scaled temperature into the relation of pressure as a function 
of energy density e by means of the above quasi-particle model provides a unique equation of state 
at large energy densities. In this respect one may arrive at the conclusion that the equation of state 
p(e) in the form pressure vs. energy density is eventually at our disposal for large values of e. As 
the calculations exhibited in the left panel of Fig. || are for \L q = we use again our quasi-particle 
model to supplement the baryon density dependence. The successful tests of the baryon density 
dependence of two-flavour QCD in [^] and the in the previous section provide confidence in the 
reliable extrapolation to \i q > 0. 

The transition region, corresponding to temperatures of the oder of the pseudo-critical temper- 
ature, T ~ T c , however, seems not yet to be settled. We construct, therefore, a family of equations 
of state by interpolating from the unique high-density QCD part to a low-density hadron resonance 
part. The resulting equation of state is used for the hydrodynamical calculation of transverse mo- 
mentum (p±) spectra and azimuthal asymmetry V2 of various hadron species for RHIC and future 
LHC energies. Keeping initial and freeze-out conditions fixed, the spectra and V2(p±) show some 
dependence on the shape of p(e) in the region around T c . This may enable further constraints on the 
equation of state by comparison with experimental data from RHIC, a programme already pursued 
by [|ll|]. For results and details the interested reader is referred to [12, |13| ]. 
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4. Including the critical point 

The change of the slope of the employed effective coupling ( |2.3| ) at T c is primarily responsible 
for the pronounced structures observed in C4, c\, C4" and cf in Fig. While the rise for T > 
T c close to T c is dictated by the perturbative behavior of G 2 , the decrease for T < T c is due to 
the conversion into the linear temperature dependence. Such a change in the curvature behavior 



was also found by solving Dyson-Schwinger equations in Coulomb gauge Q14p. This behavior 
may be interpreted as an indication of some criticality in agreement with lattice QCD results, 
nonetheless, one can further extend the model by including explicitly a singular part which belongs 
to the universality class of the three-dimensional Ising model. In such a way the conjectured QCD 
critical point can be modelled without destroying the agreement with the available lattice QCD 
results [J7]|. Such a phenomenological procedure, first proposed in Jl5|], is described in [16]. It 
offers the opportunity to study various observables within a hydrodynamical framework for matter 
states in the vicinity of the critical point, similar to first investigations along this line in Jl7|]. 

5. Recent developments 

The severe approximations described in beginning of section |2] are matter of ongoing investi- 
gations. Here we describe some of such studies. 
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Figure 2: Left panel: Scaled pressure as a function of scaled temperature (quasi-particle model [curves] 
adjusted to various selected lattice QCD data [symbols] from ||, |l(J). Right panel: Resulting equation of 
state in the form of pressure p vs. energy density e when using the scale T c = 170 MeV. The curves are robust 
against variations in T c for small and large e. Only in the transition region 1 GeV/fm 3 < e < 5 GeV/fm 3 
noticeable differences of at most 20% arise when modifying T c by ±10 MeV. 



5.1 Chiral extrapolation 

Lattice performances require still fairly large quark masses. Given the form of the approxi- 
mated quark dispersion relation (2.2) with lattice mass parameter m q one can attempt to perform a 



chiral extrapolation by putting m q — > 0. The result of such a procedure is exhibited in Fig. [| The 
semi-quantitative agreement with new lattice QCD results with smaller quark masses supports the 
idea that an implicit dependence of G 2 on quark masses is weak. 

A calculation of the quasi-particle dispersion relations based on one-loop self-energies with 



finite quark masses in Feynman gauge Jl9| ] renders the above simplified expressions ( |2.1[ ) and ( [2/2| ) 
at n q = into 

co 2 =k 2 + mi, mi= 1 -G 2 T 2 (N c + l -^(^)Y (5.1) 



!c +nr ! + 2Ml, M\ = \ G 2 T 2 (\ + ) (5.2) 



with the asymptotic representation for small values of m q /T 



a 2 = ^ (log 71 + --7 £ ), (5.4) 

T— ! 



J*{x) = 1 + ci2X 2 + aix 2 log* 2 + a$x A H , (5.3) 

3 /, 1 
^ (l0g7I+ 2 

a L = (5.5) 

«4 = ~^m. (5-6) 
The asymptotic masses m 2 and M£ = m 2 + 2M\ are found to be gauge invariant quantities [[H]] . 



(5.1) and (5.2) as well as improved dispersion relations offer the possibility of sound chiral extrap- 



olations. 
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Figure 3: Scaled pressure as a function of scaled temperature. Symbols depict lattice QCD data for 2+1 
flavors from [[| (full circles, m u j = OAT, m s = T) and jl8| ] (squares, m n = 220 MeV). The black solid 
curve employs quark masses as on the lattice, the red dashed curve is for the chiral limit, while the blue 
dash-dotted curve uses m ud = 0.0247 MeV and m s = 0.247 MeV. 



5.2 Landau damping 

While the quasi-particle model as outlined in section ^| describes fairly well the Taylor ex- 
pansion coefficients of [Q], thus establishing the correct \L q dependence for small \l q , it exhibits 
an ostensible ambiguity at larger \i q . By solving Peshier's equation to determine G 2 (T,/j, q ), the 
characteristic curves emerging from the vicinity of T c cross each other at larger pL q . As indicated 
in [p0|], including the imaginary parts of HTL/HDL self-energies cures this unpleasant feature. A 



detailed study [21] shows that the negative longitudinal gluon and (anti)plasmino contributions to 



the entropy density allow for a better adjustment of the model to lattice QCD data [18] leading to 
less crossings. However, the complete cure of crossings is due to Landau damping terms within 
Peshier's equation. 

5.3 Imaginary chemical potential 

The notorious sign problem of the fermionic determinant is avoided for a purely imaginary 
chemical potential (see [22] for recent reviews), where QCD recovers the center symmetry giving 
rise to the Roberge-Weiss periodicity [^3|]. Our quasi-particle model, as described in section || 
can be applied accordingly by the replacement of baryo-chemical potential jj, q — > = /ig/3. This 
replacement flips signs at a few important points in the equations for the thermodynamic quantities, 
the self-energies and the Peshier equation. For instance, and M\ in ( jlH] ) and (2.2) render to 



mi 



1 ( [2N C + Nf] T 2 - ^Nftf) G 2 (T, ifr), (5.7) 

This comprises an additional sensible test, even if the truncation of QCD by the chosen <I> func- 
tional, the approximated HTL/HDL self-energies and the neglect of imaginary parts in the self- 
energies discards the center symmetry of full QCD with imaginary chamical potential. A compari- 
son with lattice QCD results is displayed in Fig. ||. In particular, the importance of the BKS-effect 
[Q, ^] for the successful description of the observed pattern in the quark number density can be 
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Figure 4: Left: scaled quark number density as a function of imaginary chemical potential for various 
temperatures. Symbols denote lattice QCD data from [ p4| for T = 1.1, 1.5, 2.5, 3.5 T c (diamonds, circles, 
squares and triangles, respectively). Right: test of the BKS effect for T = 1.1 T c as explained in the text. 



studied as depicted in the rigth panel of Fig. ^ For example, one could neglect the terms explic- 



itly depending on jit,- in the quasi-particle dispersion relations with (5.7) and (5.8) which enter the 
thermodynamic expression of the quark number density but leaving Peshier's equation unchanged. 
This renders the found results only at larger values of ]X[ (dashed line). Thermodynamic consis- 
tency, nonetheless, requires in addition corresponding changes in Peshier's equation resulting in 
negligible deviations (dash-dotted line) from the results shown in the left panel of Fig. Q Thus, 
one may conclude that an explicit /x, -dependence of the quasi-particle's effective masses is not 



necessarily required for describing these lattice QCD results p4| ] ; instead the /i, -dependence of G 
matters. 



6. Strongly interacting matter in the expanding universe 

The WMAP data on the fluctuations of the cosmic microwave background deliver a value of 
6.1 x 10~ 10 (ACDM model and 3-year WMAP-data-only [[25|]) for the ratio of baryons to photons. 
Assuming adiabaticity and baryon conservation this translates into 2.3 x 10~ n for the inverse of 
the specific entropy of baryons. The corresponding adiabatic path is displayed in the temperature- 
baryo-chemical potential plane in Fig. || (for plots in log scales cf. [26,^7]). A similar path has been 



reported in [28], where also charge neutrality and lepton conservation are implemented. Intriguing 



is the sharp turn from a region of small baryo-chemical potential, /ig /T <C 1 , to large one, }Xb/T ^> 
1, at temperature scale slightly below 50 MeV. On the displayed path the time varies from 1 /isec 
to 0. 1 sec, and the universe is radiation dominated. At T > T c the strongly interacting (deconfined) 
matter dominates over the electro-weak matter, while at T < T c the electro-weak matter dominates 
for a long time until recombination. "Dominating" means here that the contribution to energy 
density, pressure or entropy density exceeds the other contributions. 



7. Summary 

In summary we survey the status of our quasi-particle model and compare it with lattice QCD 



7 



Do we know eventually p (e) ? 



B. Kampfer 



250 



200 - 



> 




P B [MeV] 



Figure 5: Adiabatic path of cosmic matter for inverse specific entropy of 2.3 x 10 , assuming about 10 
degrees of freedom (photons, standard model neutrinos, electrons) for the entropy density. The asterisks 
depict chemical freeze-out points from (table 2-upper part, and a LHC estimate mentioned in the text 
there). 

data. We argue that the equation of state in the form pressure as a function of energy density, pie), 
is eventually known at large energy densities and small baryo-chemical potential. 
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